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ABSTRACT 



An experimental investigation was carried out to determine the effect on the 
convective heat transfer coefficient of protuberances in a rectangular passage. 

A cross-flow compact type heat exchanger formed by several parallel spaced 
plates including protuberances was tested by using water, glycol and air as a working 
fluids. Heat transfer coefficient data was obtained by measuring mass flow-rate, 
temperature and pressure drops across the heat exchanger. 

Heat transfer characteristics are discussed and correlations for Nusselt number, 
Colburn factor and friction factor have been obtained. 
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NOMENCLATURE 



A : Heat transfer surface area 
A^ : Bump base area 
A^ c : Bump cross sectional area 
A^ s : Bump surface area 
Ap : Passage flow area 
A g : Plate surface area 
b : Bump base diameter 
C c : Cold fluid capacity rate 
C^ : Hot fluid capacity rate 
C r : Capacity rate ratio 
Cp : Specific heat at constant pressure 
D : Equivalent diameter 
f : Friction factor 
G : Mass velocity 
h : Heat transfer coefficient 
: Bump height 
j : Colburn factor 
k : Thermal conductivity 
K c : Contraction loss coefficient 
K e : Expansion loss coefficient 
L : Plate length 
m : Fluid mass flow rate 
N : Bump numbers 
N tu : Number of transfer units 
Nu : Nusselt number 
P : Wetted perimeter 
Pr : Prandtl number 
Q : Heat tranfer from hot fluid to cold 
q : Heat flux 
Re : Reynolds number 
r w : Wall resistance 
s : Arc length of the bump 
St : Stanton number 



[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m 2 ] 

[m] 

[J/kgs] 

[J/kgs] 

[J/kgK] 

[m] 

[kg/m 2 s] 

[W/m 2 K] 

[m] 

[W/mK] 



[m] 

[kg/s] 



[m] 

fluid [W] 

[W/m 2 ] 

[W/m 2 K] 

[m] 
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T : Hot fluid temperature 


[C] 


t : Cold fluid temperature 


[C] 


: Hot fluid bulk temperature 


[C] 


t c : Cold fluid bulk temperature 


[C] 


T : Hot fluid average temperature 


[C] 


t : Cold fluid average temperature 


[C] 


T s : Plate surface temperature 


[C] 


tp : Passage height 


[m] 


U : Overall heat transfer coefficient 


[W/nrK] 


V : Fluid velocity 


[m/s] 


v : Specific volume 


[m 3 /kg] 


w : Plate width 


[m] 


AT : Temperature drop 


[C] 


p : Fluid density 


[kg/m 3 ] 


5 : Material thickness 


[m] 


u : Kinematic viscosity 


[m 2 /s] 


p : Dynamic viscosity 


[Ns/m 2 ] 


p w : Dynamic viscosity at wall temperature 


[Ns/m 2 ] 



£ : Exchanger effectiveness 

<7 : Ratio of free flow area to frontal area 



Subscripts : 

c : coldjluid 
h : hot fluid 

1 : inlet 

2 : outlet 
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I. INTRODUCTION 



A. BACKGROUND 

In a heat exchanger, one of the important parameters controlling the net heat 
transfer from the hot fluid to the cold fluid is the surface area separating the two fluids 
across which the heat transfer takes place. When large quantities of heat are to be 
transferred, computations usually indicate a requirement for large heat transfer surface 
areas. Increasing the area necessarily means either increasing the total path length 
travelled by the fluids in the heat exchanger or simultaneously decreasing the diameter 
of tubes and increasing the number of tubes. 

Where high volume is not a limiting constraint, flat plates can be employed. 
However, use of flat plates to obtain the necessary surface may involve other heat 
transfer considerations in the device. 

In order to enhance the heat transfer characteristics of flat plate exchangers, 
some type of obstruction may be placed between the plates in order to cause the fluids 
to take a sinuous path between them. These obstructions are often in the form of a 
honeycomb structure but they inherently possess several problems. First, there is 
usually a straight path through the heat exchanger in the direction of flow thereby 
defeating the purpose of the obstructions. In addition, the sharp bends and corners in 
the honeycomb structure can create undesirable dead spots within the heat exchanger. 
They may also represent weak points and areas of deleterious stress concentration in 
the structure which, upon heating and cooling, will often crack. 

Interest in heat exchanger surfaces with a high ratio of heat transfer area to core 
volume is increasing at an accelerated pace. The primary reasons for the use of these 
more compact surfaces is that smaller, lighter-weight, and lower cost heat exchangers 
are the result. These gains are brought about both by the direct geometric advantage 
of higher 'area density' and also a higher heat transfer coefficient for the smaller flow 
passages. 

B. THESIS OBJECTIVES 

The objective of this thesis is to determine the effect of protuberances in a 
rectangular passage on the convective heat transfer coefficient. 
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In order to accomplish this objective, measurements of fluid mass flow-rate, inlet 
and outlet temperatures and exchanger pressure drops were obtained for several 
different fluids flowing in a cross-flow compact type heat exchanger containing the 
protuberances. Heat transfer coefficients were then obtained. From these data, 
correlations for Nusselt number, Colburn factor and friction factor were then 
developed. 
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II. THEORY OF OPERATION 



A. GENERAL 

The rate of heat transfer between a solid surface and a fluid can be expressed in 
much the same way as for simple heat conduction through solids: 

Q = h A AT (eqn 2.1) 

Here AT is the temperature drop between the bulk fluid and the wall, A is the heat 
transfer surface area and h is the heat transfer coefficient. The last term is equivalent 
to the term kjL in the simple heat conduction equation. 

Equation (2.1), which is often referred to as "Newton's Law of Cooling" is easy to 
use except for the evaluation of the heat transfer coefficient h, which depends on the 
fluid flow conditions, the thermal properties of the fluid, and the passage size. 

It is often convenient in heat exchanger work to combine the two surface heat 
transfer coefficients (that is, for the hot and cold fluid streams) and the thermal 
resistance of the tube wall to give a single parameter: the overall heat transfer 
coefficient, U, which may be defined as 

U = Q / A AT (eqn 2.2) 

B. INDIVIDUAL AND OVERALL HEAT TRANSFER COEFFICIENTS 
Equation (2.2) may be written in another form : 

q" = U (Tj - T 2 ) (eqn 2.3) 

Here q is the heat flux from stream- 1 to stream-2, having bulk temperatures T [ and T 2 
and U is the overall heat transfer coefficient. 

It is now time to consider the questions of what precisely is meant by the 
interface area, and what is the relation of U to the individual resistances to heat 
transfer represented by the solid wall separating the two fluids, the deposits of dirt that 
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may cling to the surfaces of the wall, and the low-velocity fluid films that adjoin those 
deposits or the wall itself. 

Fig. 2.1 illustrates the temperature profile discontinuities at locations a (inner 
surface of inner deposit), b (outer surface of inner deposit and inner surface of tube 
wall), c (outer surface of tube wall and inner surface of outer deposit), d (outer surface 
of outer deposit). It can also be seen that the fluid temperatures lie both above and 
below T a and 7^, for these are defined as flow averages, that is, averages for the duct 
cross sections in which each element of area is weighted by the specific heat C p times 
the local mass flow rate pV. 

Thus, 

T lor2 = ( f P C p V T 3s / J p C p V ds ) 10r2 (eqn 2.4) 

where the averages are taken over the whole of the cross section of the ducts in 
question. 




Figure 2.1 Temperature distribution in the neighborhood of a tube wall. 

Let h^ a and now represent the film heat transfer coefficient on the inside and 
outside of the tube, respectively; let h and h cc f represent heat transfer coefficients for 
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the two dirt deposits; and hf JC represents the heat transfer coefficient for the tube wall 
itself. h a fo, hfo c , and h will all be related to the relevant material thickness 5 and 
thermal conductivities k by 



^ab ^ab ' ^ab 



(eqn 2.5) 



^bc ^bc 1 ^bc 



(eqn 2.6) 



^cd ^cd ' ^cd 



(eqn 2.7) 



Bach heat transfer coefficient, to be meaningful, must be associated with an area 
A. When this is done, the fact that the heat flow through all the elements in question 
is the same leads to 



Here A Q ^, A^ c , A C( j can most appropriately be taken as the arithmetic means of 
the inner and outer surface areas of the materials in question; and A [a and A ^ arc 
best taken as the surface areas pertaining to locations a and d respectively. The area A 
without subscript, associated with U, requires arbitrary definition; it can be set equal to 
any of the other areas. 

Elimination T a , 7^, T , and T ^ now leads to the following important connection 
between the overall coefficient U and the individual coefficients 



UA(T, -T 2 ) = h la A |a (T, -T a ) 



(eqn 2.8) 



- ^ab A ab ^a ~~ ^b^ 

- ^bc A bc ^b ~~ T c) 
= h cd A cd ( T c “ T d^ 
= ^d2 A d2 ^d — ^ 2 ) 



U - [( A / h Ja A la ) + ( A / h ab A ab ) + ( A / h bc A bc ) 
+ ( A / h cd A cd ) + ( A / h d2 A d2 )] “ 1 



(eqn 2.9) 
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Evidently, the relationship is best understood as implying that the reciprocal of 
UA is the overall resistance to heat transfer, and that the overall resistance is equal to 
the sum of the individual resistances (/t la /l la ) ^A^ab^al) anc ^ 

In the case of heat transfer across a plane wall, all the areas take identical values. 
Then Equation (2.9) reduces to 

U = [( 1 / h la ) + ( i / h ab ) + ( 1 / h bc ) + ( I / h d2 )]“ 1 (eqn 2.10) 

C. HEAT TRANSFER COEFFICIENT DEPENDENCES 

A dimensional analysis easily shows that the heat transfer coefficient may be 
found to be related to two dimensionless parameters : the Nusselt and Prandtl 
numbers. The Nusselt number, Nu, is the term used for the grouping hLjk, where L is a 
typical length dimension that pertains to the apparatus, and k is the thermal 
conductivity of the fluid. This quantity is the ratio of the heat transfer coefficient to the 
thermal conductance implied by the quantity kj L. It can be seen intuitively, from the 
simple conduction equation, that the heat flow rate to a fluid flowing through a 
passage should be proportional to the thermal conductivity divided by a representative 
distance in the direction of heat flow such as the passage equivalent diameter. 

Obviously, further prescription is needed of precisely what dimension is in 
question (such as duct diameter), and just where the thermal conductivity is to be 
measured (such as at the interface, at the bulk condition, or in some reference state). 

The Prandtl number, Pr, is the term used for the ratio C^\ t/k This quantity is the 
ratio of the molecular diffusivity of momentum (as indicated by the viscosity) to the 
molecular diffusivity of heat (as indicated by the ratio of the thermal conductivity to 
the specific heat). 

D. HEAT TRANSFER COEFFICIENT LIMITATIONS 

The concept of the heat transfer coefficient is so convenient, and so reminiscent 
of electrical circuit theory, that there is some danger that it's validity may be presumed 
to be unlimited. A few remarks on it's limitations are therefore in order. 

First, the heat flux (to or from the fluids) may not be proportional to the 
temperature difference alone. Other examples can be found, for example, those 
involving chemical reaction rather than phase change. 
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Second, the value of the coefficient is by no mean independent of the values of 
the property differences that multiply it. In the case of heat transfer by natural 
convection, it is well known that the coefficient increases with the imposed 
temperature difference. 

Third, there are circumstances in which nonuniformities of temperature or 
concentration along the interface are so steep as to affect the value of the transfer 
coefficient. In these circumstances, it may be best to work with the fluxes and 
property differences directly, avoiding all reference to coefficients. 
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III. THE EMPRICIAL CORRELATION FOR THE CONVECTIVE HEAT 

TRANSFER COEFFICIENT 



The manner in which a convection heat transfer coefficient may be obtained 
experimentally is illustrated in Fig. 3.1. In a prescribed geometry, such as a flat plate in 
parallel flow, convection heat transfer occurs from the plate to the fluid, if the surface 
of the plate is at a higher temperature than that of the fluid. 




Figure 3.1 Experiment for measuring the average convection 

heat transfer coefficient. 

The convection heat transfer coefficient h, which is an average associated with 
the entire plate, can be computed from Newton's law of cooling, 

Q=hA s (T $ -T b ) (eqn 3.1) 

where 



h : average convective coefficient 


[W/m 2 K] 


A $ : plate surface area 


[m 2 ] 


T $ : plate surface temperature 


[C] 


T^ : fluid bulk temperature 


[C] 
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From a knowledge of the characteristic length L and the fluid properties, the Nusselt, 
Reynolds, and Prandtl numbers can be computed from their definitions. The Nusselt 
number is a function of the Reynolds and Prandtl numbers. 

Nu = ( h L / k ) = f( Re, Pr ) (eqn 3.2) 

and 

Re = VL/u . (eqn 3.3) 



Pr = Cp p / k (eqn 3.4) 

where 

L : characteristic length [m] 

k : thermal conductivity [W/mK] 

V : fluid velocity [m/s] 

p : fluid dynamic viscosity [Ns/m 2 ] 

Cp : specific heat [J/kgK] 

0 : fluid kinematic viscosity [m 2 /s] 



The forgoing procedure is reproducible for a variety of test conditions. One can vary 
the velocity V and the plate length L, as well as the nature of the fluid, using for 
example air, water and glycol, which have substantially different Prandtl numbers. We 
would then be left with many different values of the Nusselt numbers corresponding to 
a wide range of Reynolds and Prandtl numbers and the results could be plotted on 
log-log scale, as shown in Fig. 3.2. 

Each symbol represents a unique set of test conditions. As is often case, the 
results associated with a given fluid, and hence a fixed Prandtl number, fall close to a 
straight line, which may be represented by an algebraic expression of the form 
Since the values of C, m and n are independent of the nature of the fluid, the family of 
straight lines corresponding to different Prandtl numbers can be collapsed to a single 
line by plotting the results in terms of the ratio, Nu/Pr n , as shown in Fig. 3.3. 
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Figure 3.2 Dimensionless representation of convection 
heat transfer measurements. 



Nu = C Re m Pr 11 (eqn 3.5) 

Because Equation (3.5) is inferred from experimental measurements, it is termed an 
empricial correlation. However, the specific values of coefficient C and the exponents m 
and n vary with the nature of the surface geometry and the type of flow. 
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Figure 3.3 Dimensionless representation of convection 
heat transfer measurements. 
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IV. EXPERIMENTAL PROCEDURE 



A. EXPERIMENTAL PROCEDURE 

To reach the thesis objectives an experimental system was established. It is 
shown schematically in Fig. 4.1, and using water, glycol and air as a working fluids, the 
effect upon the convective heat transfer coefficient by the protuberances was 
investigated. 

The physical configuration of the heat exchanger and some instrumentation 
methods will be discussed in the following sections. 

1. Heat Exchanger 

The heat exchanger (Fig. 4.2) was formed by 14 parallel spaced plates with the 
spaces between the plates defining fluid passages. The plates are totally nested with 
each other to provide, by minimizing unsupported areas, a structurally sounder device. 

When the fluids enter the passage, they are confronted with a unique pattern 
of interfering protuberances which come downward from the plate above and upward 
from the plate below. The turbulence created as the fluids make their sinuous path 
from inlet to outlet greatly enhances the heat transfer characteristics. Moreover, the 
fact that the rounded protuberances nest within each other and actually engage the 
adjacent plate gives the heat exchanger structural rigidity which permits the use of 
thinner metallic plate material. This further enhances the heat transfer characteristics 
by lowering the metal resistance to heat flow. 

A brief description of the figures now follows : 

Fig. 4. 2. a is an external view of the heat exchanger. 

Fig. 4.2.b is an internal view of the heat exchanger. 

Fig. 4.3 is a detached edge view of the top plate 

configuration employed in the heat exchanger. 

Fig. 4.4, and Fig. 4.5 arc detached edge views of another 

type of configuration employed in the heat exchanger. 

Fig. 4.6 is a detached edge view of the bottom plate 
configuration employed in the heat exchanger. 

Fig. 4.7 is a view taken substantially along line a-a of 
Fig. 4.2, and shows the configurations of the plates 
in the heat exchanger. 
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Figure 4. 1 Schematic diagram of the experimental system. 
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2. Geometric Properties Of The Heat Exchanger 



PLATE 

Plate surface area : A $ = Lw 

Bump surface area ’ ^bs 2 it r 

Bump cross section area : A^ c = (7t r~ <I>/360) - (bd/2) 

Bump base area : A^ = K b 2 /4 




where 



A $ : plate surface area 


[m 2 ] 


L : plate length 


[m] 


w : plate width 


[m] 


A| 3S : bump surface area 


[m 2 ] 


A^ c : bump cross section area 


[m 2 ] 


A^ : bump base area 


[m 2 j 


h^ : bump height 


[m] 


b : bump base diameter 


[m] 
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Figure 4.2 Schematic view of the heat exchanger. 
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Figure 4.4 A detached edge view of another type of configuration. 





Figure 4.6 A detached edge view of the bottom plate configuration. 
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Figure 4.7 The plate configurations in the heat exchanger. 

PASSAGE FLOW AREAS 

Hot passage flow area : A^ = L t p - A 

Cold passage flow area : Ap c = w t - N c A 

where 

tp : passage height, [m] 

and N' c are the numbers of bumps on hot and cold sides of the plate, 

respectively. 

EQUIVALENT DIAMETER 
D = 4A f /P 

where 

Ap: passage flow area 
P : wetted perimeter 

P h = 2 L + N h ( s - b ) 

P c = 2 w + N c ( s - b ) 

where 

s : arc length of the bump 
b : bump base diameter 
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[m 2 ] 

[m] 



[m] 

[m] 



: hot wetted perimeter 
P c : cold wetted perimeter 



[m] 

[m] 



3. Instrumentation 

Temperature measurements were made using Coppcr-Constantant 
thermocouples using a distilled water-ice bath reference. A total of 16 thermocouples 
were used with individual selection being made with a thermocouple selector switch. 
This switch and an individual connector plug were installed iri a tightly closed, 
insulated metal box to keep all the connecting junctions at a uniform temperature and 
thus reduce errors from connecting junctions. Four thermocouples were used on each 
inlet and outlet section of the heat exchanger. 

Pressure drops across the heat exchanger were measured using tw r o mercury 
manometers. 

Other instrumentation included digital flow-meters used to measure and 
control mass flow-rate so that a steady-state condition was assured. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The effects on the convective heat transfer coefficient of protuberances has been 
investigated in the fully-developed, laminar flow regime. 

The effects of protuberances are beneficial as they : 

1. Increase the structural rigidity which permits the use of thinner plate material 
and further enhances the heat" transfer characteristics. 

2. Create turbulence which augments the heat transfer characteristics. 

The heat transfer correlations obtained are as follows : 

Nu = 0.113 Re 62 Pr 1 / 3 (H/H w )- 14 

j = 0.109 Re"' 38 0i/M w )- 14 

f = 0.549 Re~ ,4 ° (|x/|J. w ) ~ ,14 

B. RECOMMENDATIONS 

The present study has been associated with the laminar llow regime. Additional 
testing should be carried out in the turbulent flow regime. 

A data acquisition system should be employed in future testing to improve 
accuracy. 
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APPENDIX A 

SAMPLE CALCULATIONS 



CONSTANT VALUES 

Heat exchanger overall dimensions : 
length 
width 
height 
Plate length 
Plate width 
Plate area 
Bump surface area 
Bump cross sectional area 
Bump base area 
Effective heat transfer area 
(per plate) 

Total effective heat 
transfer area (one side) 

Hot passage flow area 

Cold passage flow area 

Hot side equivalent diameter 

Cold side equivalent diameter 

Flow passage height 

Bump arc length 

Bump arc radius 

Bump base diameter 

Hot wetted perimeter 

Cold wetted perimeter 

Bump numbers on the hot side 

Bump numbers on the cold side 



209.55mm 


(8.25in) 


205.58mm 


(8.094in) 


63.5mm 


(2.5in) 


108.74mm 


(4.28125in) 


104.78mm 


(4.125in) 


1 1393.8mm 2 


(17.66in 2 ) 


108.7mm 2 


(0.16855in 2 ) 


226mm 2 


(0.3503in 2 ) 


77.1mm 2 


(0. 1 195in 2 ) 


14243.9mm 2 


(27.078in 2 ) 


199414.2mm 2 


(309.092in 2 ) 


141.9mm 2 


(0.21989in 2 ) 


106.7mm 2 


(0.1653in 2 ) 


2.36mm. 


(0.093in) 


1.82mm. 


(0.01 \ 55m) 


3. 175mm. 


(0.125in) 


12.426mm. 


(0.489in) 


5.45mm. 


(0.215in) 


9.906mm. 


(0.39in) 


240.1S2mm. 


(9.456in) 


234.752mm. 


(9.242in) 



9 

10 
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2. SAMPLE CALCULATIONS 

In this thesis, the following calculation procedure was used to calculate heat 
transfer characteristics and to obtain the correlations. 

a. CALCULATION PROCEDURE 

1. Determine the bulk temperatures, T^ and t c , 

T b = ( T i + T 2 > / 2 

and 

t c = ( + t 2 ) / 2 

2. Determine p, Cp, p and k and the Capacity rates, C c and C^, 

C c - m c^pc 

and 

C h = m h C ph 

3. Determine the exchanger effectiveness, c, 

£= C h (T, - T 2 )/ Cmin (Ti - t,) ifC h <C c 
or 

c - c c ( t 2 “ tj ) / Tj — tj ) if C c <C h 

4. Determine the number of transfer units, N tu , from C-N tu relations, both fluids 
unmixed, ( Ref.6. ) 

c = 1 - exp[(l/C r )(N lu )' 22 [cxp(- C r (N tu )' 78 )- 1 ]] 
where C f is the capacity rate ratio, C ni j n /C max . 

5. The overall heat transfer coefficient, U, is formed from the number of transfer units. 

^ “ ^tu^min ! ^ 

where A is the total effective heat transfer area. 

6. The convective heat transfer coefficient, h^, will be 

h h = Id/U) - (l/h c ) - r w ] _1 

Where h^ is the hot fluid heat transfer coefficient 
h c is the cold fluid heat transfer coefficient 
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r,., is the wall resistance 

7. The fluid velocity, V, is for either fluid 

V = m / pAf 

where Apis the passage flow area. 

8. The wall temperature, T w , is computed from 

T w =7+ [h h /(h h + h c )](T -T) 

where T is the average temperature of the hot fluid, 
t is the average temperature of the cold fluid. 

10. Form the Prandtl number, 

C p p/k 

the Reynolds number, 
pVD/n 

and the viscosity correction factor, 

with p w calculated at wall temperature. 

1 1. The friction factor, f, is determined from 

f = (A c v m /Av 1 )l(2g c AT/G J v 1 ) + (l- K^Vj/V;)- (K c +1- <T 2 )- 2(vjy [ — 1)] 

where v m : fluid specific volume at mean condition. 

Vj : fluid specific volume at inlet condition. 
v 0 : fluid specific volume at outlet condition. 

G : mass velocity. 

A : total heat transfer area. 

Ap : free flow area. 

<7 : ratio of free-flow area to frontal area. 

<7^ = 0.028 for hot side. (Ref. 6.) 

<7 C = 0.0178 for cold side. (Ref. 6.) 

K c : contraction loss coefficient. 
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K c = 1.3 both hot and cold side. (Ref. 6.) 

K p : expansion loss coefficient. 

K e k = 0.92 for hot side. (Ref. 6.) 

K £C = 0.97 for cold side. (Ref. 6.) 

b. HOT WATER-COLD WATER STUDY 
The experimental data is tabulated in Table 1. 

The experiment results are tabulated in Tables 2, 3 and 4. 

c. HOT GLYCOL-COLD WATER STUDY 
The experimental data is tabulated in Table 5. 

The experimental results are tabulated in Tables 6, 7 and 8. 

d. HOT WATER-COLD AIR STUDY 

The experimental data is tabulated in Table 9. 

The experimental results are tabulated in Tables 10 and 11. 

e. CALCULATION OF CORRELATIONS 

1. THE CORRELATION FOR THE NUSSELT NUMBER : 

All of the experimental data for a correlation for the Nusselt number are 
plotted in Figure A.l. Observe that the results show very little scatter, because of this, 
one may take runs 5 and 8 in the water-water study and begin with a proposed form : 

Nu = C Re m Pr n (n/n w )’ 14 

Then, it is only a matter of algebra to show that 

LogNu = LogC + mLogRe + nLogPr + 0. 14Log(n/n w ) 

and from the Table 4, for runs 5 and 8 respectively 

Log(10.2) = LogC + mLog(63 1.8) + (l/3)Log(4.79) + 0.14Log(0.S92) 

Log(l 1.4) = LogC + mLog(756.4) + (l/3)Log(4.S0) + 0. 14Log(0.887) 

from which one may solve for C and m, 

C = 0.113 
m = 0.62 
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Hence, the correlation is 

Nu = 0.113 Re’ 62 Pr 1 ' 3 (p/p w )' 14 

2. THE CORRELATION FOR THE COLBURN FACTOR : 

All of the experimental data for a correlation for the Colburn factor are 
plotted in Figure A. 2. Observe that the results show very little scatter. Because of this, 
one may take runs 5 and 15 in the glycol- water study and begin with a proposed form : 

j = C Re m (n/n w )' 14 and 

Logj = LogC + mLogRe + 0. 14Log(p/p w ) 

and from the Table 8, for runs 5 and 15 respectively 

Log(0.025S) = LogC + mLog(42.77) + 0.14Log(0.651) 

Log(0.0200) = LogC + mLog(83.15) + 0.14Log(0.652) 

from which one may solve for C and m, 

C = 0.109 
m = — 0.38 

Hence, the correlation is 

j = St Pr 2 / 3 (H/[l , v )- 14 

j = 0.109 Re - ' 38 (|i/n w )- 14 

where St is the Stanton number, St = h/pCpV. 

3. THE CORRELATION FOR THE FRICTION FACTOR : 

All of the experimental data for a correlation for the friction factor arc plotted 
in Figure A. 3. Observe that the results show very little scatter. Because of this, one 
may take runs 2 and 18 in the water- water study and begin with a proposed form : 

f = C Re m 0i/p w )“- 14 and 

Logf = LogC + mLogRe - 0. 14Log(p/p w ) 

and from the Table 11, for runs 2 and 18 respectively 
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Log(0.051) = LogC + mLog(378.4) - 0.14Log(0.896) 
Log(0.03l) = LogC + mLog( 1301.0) - 0.14Log(0.897) 

from which one may solve for C and m, 

C = 0.549 
m = - 0.40 

Hence, the correlation is 

f = 0.549 Re - ' 40 (n/n w )“- 14 
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TABLE 1 

THE EXPERIMENTAL DATA FOR THE WATER- WATER STUDY 



RUN 


AP 


m h 


rru 


HOT WATER 


COLD WATER 


NO 


inH 2 0 


c 


T l 


T 2 




4 


1 


1.36 


0.089 


0.0S9 


37.38 


33.06 


21.90 


26.22 


2 


2.04 


0.114 


0.114 


37.40 


33.42 


22.82 


26.80 


3 


3.13 


0.146 


0.146 


37.30 


33.71 


22.94 


26.53 


4 


4.08 


0.171 


0.171 


37.28 


33.72 


22.87 


26.43 


5 


4.76 


0.190 


0.190 


37.25 


33.73 


22.70 


26.22 


6 


5.44 


0.206 


0.206 


37.15 


33.66 


22.55 


26.04 


7 


5.71 


0.216 


0.216 


37.20 


33.75 


22.39 


25.84 


8 


6.12 


0.229 


0.229 


37.10 


33.67 


22.12 


25.55 


9 


7.48 


0.248 


0.248 


36.90 


33.52 


21.88 


25.26 


10 


8.16 


0.267 


0.267 


36.80 


33.48 


21.85 


25.17 


11 


9.65 


0.2S6 


0.286 


36.10 


32.96 


21.70 


24.84 


12 


10.47 


0.305 


0.305 


35.90 


32.77 


21.41 


24.54 


13 


12.24 


0.330 


0.330 


35.50 


32.44 


21.23 


24.29 


14 


13.32 


0.349 


0.349 


35.30 


32.30 


21.15 


24.15 


15 


13.60 


0.368 


0.368 


35.00 


32.13 


21.18 


24.06 


16 


15.64 


0.381 


0.381 


34.60 


31.80 


20.81 


23.61 


17 


16.32 


0.400 


0.400 


34.50 


31.75 


20.88 


23.63 


18 


17.68 


0.412 


0.412 


34.35 


31.68 


21.00 


23.67 
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TABLE 2 










PROPERTIES OF THE HOT WATER 




RUN 

NO 


T bh 


P 


S 


\i 10 6 


k 10 3 


1 


35.22 


1006.30 


4178 


721.344 


625.152 


2 


35.41 


1006.40 


4178 


718.532 


625.456 


J 


35.51 


1006.40 


4178 


717.126 


625.608 


4 


35.50 


1006.40 


4178 


717.200 


625.600 


5 


35.49 


1006.40 


4178 


717.348 


625.584 


6 


35.41 


1006.36 


4178 


718.606 


625.448 


7 


35.48 


1006.40 


4178 


717.570 


625.560 


8 


35.39 * 


1006.35 


4178 


718.902 


625.416 


9 


35.21 


1006.30 


4178 


721.492 


625.136 


10 


35.14 


1006.26 


4178 


722.528 


625.024 


11 


34.53 


1006.01 


4178 


731.556 


624.048 


12 


34.34 


1005.90 


4178 


734.442 


623.736 


13 


33.97 


1005.80 


4178 


739.844 


623.152 


14 


33.80 


1005.70 


4178 


742.360 


622.880 


15 


33.56 


1005.60 


4178 


745.875 


622.500 


16 


33.20 


1005.50 


4178 


751.240 


621.920 


17 


33.13 


1005.40 


4178 


752.350 


621.800 


18 


33.02 


1005.40 


4178 


753.980 


621.624 
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TABLE 3 

SOME RESULTS OF THE WATER-WATER STUDY 



RUN 

NO 


£ 


*tu 


U 


h h 




T 

w 


1 


0.279 


0.396 


737.40 


1510.2 


1510.2 


29.64 


2 


0.273 


0.384 


919.48 


1894.2 


1894.2 


30.11 


3 


0.250 


0.339 


1037.25 


2145.0 


2145.0 


30.12 


4 


0.247 


0.334 


1199.68 


2494.2 


2494.2 


30.08 


5 


0.242 


0.324 


1293.06 


2696.7 


2696.7 


29.98 


6 


0.239 


0.319 


1379.15 


2884.4 


2884.4 


29.85 


7 


0.233 


0.308 


1393.09 


2914.9 


2914.9 


29.80 


8* 


0.229 


0.301 


1441.46 


3021.0 


3021.0 


29.61 


9 


0.225 


0.294 


1525.16 


3205.3 


3205.3 


29.39 


10 


0.222 


0.289 


1613.16 


3400.2 


3400.2 


29.32 


11 


0.218 


0.282 


1687.52 


3565.8 


3565.8 


28.90 


12 


0.216 


0.279 


17S0.73 


3774.5 


3774.5 


28.50 


13 


0.214 


0.276 


1903.12 


4050.6 


4050.6 


28.36 


14 


0.212 


0.272 


1988.80 


4245.3 


4245.3 


28.23 


15 


0.208 


0.266 


2050.83 


4386.9 


4386.9 


28.10 


16 


0.203 


0.257 


2049.47 


4425.1 


4425.1 


27.70 


17 


0.202 


0.256 


2143.50 


4599.6 


4599.6 


27.69 


18 


0.200 


0.253 


2185.54 


4696.6 


4696.6 


27.68 
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N< 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 



TABLE 4 

SOME RESULTS OF THE WATER- WATER STUDY 



Nu 


Pr 


St 


Re 


P/P w 


f 


j 


5.7 


4.82 


0.0040 


293.3 


0.891 


0.055 


0.0115 


7.2 


4.80 


0.0039 


378.4 


0.896 


0.051 


0.0112 


8.1 


4.79 


0.0035 


484.5 


0.896 


0.049 


0.0099 


9.4 


4.79 


0.0035 


568.7 


0.894 


0.048 


0.0098 


10.2 


4.79 


0.0034 


631.8 


0.892 


0.043 


0.0095 


10.9 


4.80 


0.0033 


683.2 


0.892 


0.042 


0.0094 


11.0 


4.79 


0.0032 


715.7 


0.889 


0.041 


0.0091 


11.4 


4.80 


0.0031 


756.4 


0.887 


0.040 


0.0089 


12.1 


4.82 


0.0031 


816.5 


0.886 


0.040 


0.0088 


12.9 


4.90 


0.0030 


878.0 


0.886 


0.039 


0.0087 


13.5 


4.92 


0.0030 


928.9 


0.890 


0.040 


0.0086 


14.3 


4.96 


0.0029 


986.8 


0.889 


0.038 


0.0085 


15.4 


4.98 


0.0029 


1061:1 


0.890 


0.038 


0.0085 


16.1 


5.01 


0.0029 


1118.3 


0.891 


0.037 


0.0084 


16.7 


5.05 


0.0028 


1173.8 


0.892 


0.036 


• 0.0083 


16.8 


5.06 


0.0028 


1205.4 


0.891 


0.034 


0.0081 


17.5 


5.07 


0.0027 


1262.9 


0.894 


0.033 


0.0081 


17.9 


5.12 


0.0027 


1301.0 


0.897 


0.031 


0.0080 
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TABLE 5 

THE EXPERIMENTAL DATA FOR THE GLYCOL- WATER STUDY 



RUN 


AP 


m h 


m„ 


HOT GLYCOL 


COLD WATER 


NO 


inl-Ig 


n 


c 


T l 


T 




h 


1 


0.3 


0.098 


0.089 


36.5 


33.73 


21.0 


22.8 


2 


0.4 


0.125 


0.114 


36.5 


33.80 


21.0 


22.7 


3 


0.5 


0.160 


0.146 


36.5 


33.97 


21.0 


22.6 


4 


0.6 


0.188 


0.171 


36.4 


33.93 


21.0 


22.6 


5 


0.8 


0.209 


0.190 


36.4 


33.99 


21.0 


22.6 


6 


0.9 


0.227 


0.206 


36.4 


34.03 


20.9 


22.4 


7 


1.0 


0.237 


0.216 


36.4 


34.09 


20.9 


22.4 


8 


1.1 


0.251 


0.229 


36.5 


34.24 


20.9 


22.4 


9 


1.3 


0.272 


0.248 


36.5 


34.30 


21.0 


22.4 


10 


1.4 


0.293 


0.267 


36.5 


34.33 


21.0 


22.4 


11 


1.6 


0.314 


0.286 


36.5 


34.37 


21.0 


22.4 


12 


1.8 


0.335 


0.305 


36.5 


34.44 


21.0 


22.3 


13 


2.1 


0.363 


0.330 


36.5 


34.47 


21.0 


22.3 


14 


2.3 


0.383 


0.349 


36.5 


34.50 


21.0 


22.3 


15 


2.6 


0.404 


0.368 


36.5 


34.53 ' 


21.0 


22.3 


16 


2.7 


0.418 


0.381 


36.4 


34.46 


20.9 


22.1 


17 


2.9 


0.439 


0.400 


36.4 


34.50 


20.9 


22.1 


18 


3.1 


0.453 


0.412 


36.4 


34.53 


20.9 


22.1 
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TABLE 6 

PROPERTIES OF THE HOT GLYCOL 



RUN 

NO 


T bh 


P 


c p 


H 10 2 


k 10 3 


1 


35.12 


1109.94 


2452 


1.164 


254.5 


2 


35.15 


1104.91 


2452 


1.163 


254.5 


3 


35.24 


1104.96 


2452 


1.158 


254.5 


4 


35.17 


1104.94 


2452 


1.162 


254.5 


5 


35.20 


1 104.92 


2452 


1.160 


254.5 


6 


35.22 


1105.05 


2452 


1.159 


254.5 


7 


35.25 


1 105.02 


2452 


1.158 


254.5 


8 


35.37 


1105.01 


2452 


1.152 


254.5 


9 


35.40 


1105.03 


2452 


1.150 


254.5 


10 


35.42 


1105.05 


. 2452 


1.149 


254.5 


11 


35.44 


1105.02 


2452 


1.148 


254.5 


12 


35.47 


1104.96 


2452 


1.147 


254.5 


13 


35.49 


1105.03 


2452 


1.146 


254.5 


14 


35.50 


1105.05 


2452 


1.145 


254.5 


15 


35.52 


1105.05 


2452 


1.144 


254.5 


16 


35.43 


1104.96 


2452 


1.149 


254.5 


17 


35.45 


1105.00 


2452 


1.148 


254.5 


18 


35.47 


1104.96 


2452 


1.147 


254.5 
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TABLE 7 

SOME RESULTS OF THE GLYCOL- WATER STUDY 



RUN 

NO 


£ 


N tu 


U 


h h 




T w 


1 


0.179 


0.210 


270.01 


332.30 


1510.2 


24.77 


2 


0.174 


0.204 


314.73 


382.02 


1894.2 


24.74 


3 


0.163 


0.189 


372.64 


457.53 


2145.0 


25.00 


4 


0.160 


0.185 


428.15 


525.49 


2494.2 


24.95 


5 


0.156 


0.180 


462.88 


568.88 


2696.7 


24.86 


6 


0.153 


0.176 


490.38 


602.11 


2884.4 


24.83 


7 


0.149 


0.170 


495.55 


608.56 


2914.9 


24.98 


8 


0.145 


0.165 


509.23 


624.59 


3021.0 


24.99 


9 


0.142 


0.161 


538.30 


660.49 


3205.3 


24.99 


10 


0.140 


0.159 


572.51 


703.79 


3400.2 


24.89 


11 


0.137 


0.155 


597.99 


735.23 • 


3565.8 


24.86 


12 


0.133 


0.150 


617.24 


755.59 


3774.5 


24.81 


13 


0.131 


0.147 


655.32 


801.67 


4050.6 


24.66 


14 


0.129 


0.145 . 


683.71 


836.57 


4245.3 


24.65 


15 


0.127 


0.142 


706.09 


864.60 


4386.9 


24.72 


16 


0.125 


0.140 


720.08 


884.13 


4425.1 


24.74 


17 


0.123 


0.137 


739.89 


907.08 


4599.6 


24.60 


18 


0.121 


0.135 


752.25 


921.90 


4696.6 


24.55 
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TABLE 8 

SOME RESULTS OF THE GLYCOL- WATER STUDY 



RUN 

NO 


Nu 


Pr 


St 


Re 




f 


j 


1 


3.1 


112.18 


0.00130 


21.37 


0.650 


0.142 


0.0302 


2 


4.0 


112.01 


0.00125 


25.69 


0.651 


0.115 


0.0290 


3 


4.2 


111.89 


0.00116 


32.85 


0.653 


0.090 


0.0270 


4 


4.9 


112.25 


0.00114 


38.45 


0.654 


0.090 


0.0266 


5 


5.3 


112.11 


0.00111 


42.77 


0.651 


0.094 


0.0258 


6 


5.6 


111.54 


0.00108 


46.57 


0.651 


0.090 


0.0249 


7 


5.7 


111.37 


0.00104 


48.80 


0.655 


0.090 


0.0241 


8 


5.8 


111.68 


0.00101 


51.52 


0.653 


0.090 


0.0234 


9 


6.1 


111.58 


0.00098 


55.86 


0.653 


0.087 


0.0228 


10 


6.5 


111.52 


0.00098 


60.19 


0.656 


0.084 


0.0226 


11 


6.8 


111.42 


0.00095 


64.55 


0.656 


0.082 


0.0220 


12 


7.0 


111.68 


0.00092 


68.70 


0.655 


0.081 


0.0213 


13 


7.4 


111.61 


0.00089 


74.47 


0.649 


0.080 


0.0207 


14 


7.8 


111.54 


0.00088 


78.82 


‘ 0.650 


0.080 


0.0205 


15 


8.0 


' 1 1 1.49 


0.00087 


83.15 


0.652 


0.080 


0.0200 


16 


8.2 


111.89 


0.00086 


85.71 


0.656 


0.079 


0.0200 


17 


8.4 


111.82 


0.00084 


90.05 


0.647 


0.077 


0.0195 


18 


8.6 


111.75 


0.00083 


92.97 


0.647 


0.077 


0.0193 
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TABLE 9 

THE EXPERIMENTAL DATA FOR THE WATER-AIR STUDY 



RUN 


AP 


m h 


m r 10 3 


HOT WATER 


COLD AIR 


NO 


inH 2 0 


c 


T 1 


T 


l l 


4 


1 


1.2 


0.190 


1.91 


36.8 


36.8 


19.6 


36.5 


2 


2.5 


0.190 


2.87 


36.8 


36.7 


19.6 


36.3 


3 


4.0 


0.190 


3.83 


36.8 


36.7 


19.6 


36.0 


4 


6.0 


0.190 


4.79 


36.9 


36.8 


19.6 


36.9 


5 


8.0 


0.190 


5.75 


36.9 


36.8 


19.6 


35.8 


6 


10.8 


0.190 


6.70 


36.9 


36.8 


19.8 


35.6 


7 


13.7 


0.190 


7.66 


36.9 


36.7 


19.8 


35.5 


8 


17.0 


0.190 


8.62 


36.9 


36.7 


19.8 


35.3 


9 


19.3 


0.190 


9.58 


37.0 


36.8 


19.9 


35.3 
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TABLE 10 

SOME RESULTS OF THE WATER-AIR STUDY 



RUN 

NO 


e 


N tu 


U 


h h 


h c 


1 


0.983 


4.12 


39.74 


2696.7 


40.58 


2 


0.971 


3.58 


51.88 


2696.7 


52.99 




0.953 


3.09 


59.76 


2696.7 


61.23 


4 


0.942 


2.88 


69.66 


2696.7 


71.67 


5 


0.936 


2.78 


80.72 


2696.7 


83.43 


6 


0.924 


2.61 


88.31 


2696.7 


91.56 


7 


0.918 


2.53 


97.86 


2696.7 


101.88 


8 


0.906 


2.40 


104.47 


2696.7 


109.06 


9 


0.901 


2.35 


113.69 


2696.7 


119.14 
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TABLE 11 

SOME RESULTS OF THE WATER-AIR STUDY 



RUN 

NO 


Nu 


Pr ' 


St 


Re 


f 


j 


1 


2.8 


0.707 


0.0156 


251.8 


0.055 


0.0124 


2 


3.7 


0.707 


0.0136 


378.4 


0.051 


0.0108 


3 


4.2 


0.707 


0.0118 


504.9 


0.046 


0.0094 


4 


5.0 


0.707 


0.0111 


' 631.5 


0.044 


0.0088 


5 


5.8 


0.707 


0.0108 


758.0 


0.040 


0.0086 


6 


6.3 


0.707 


0.0101 


883.3 


0.040 


0.0080 


7 


7.0 


0.707 


0.0098 


1009.8 


0.039 


0.0078 


8 


7.5 


0.707 


0.0093 


1136.4 


0.038 


0.0074 


9 


8.2 


0.707 


0.0092 


1263.0 


0.035 


0.0073 
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Figure A.l A correlation for the Nusselt number. 
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Figure A.2 A correlation for the Colburn factor. 
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Figure A. 3 A correlation for the friction factor. 
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APPENDIX B 

DIMENSIONAL ANALYSIS 



The dependent variable for the process is the convection heat transfer coefficient, 
h. For incompressible low-speed flow the independent variables that determine the 
heat transfer coefficient are the fluid velocity, V, a linear dimension (e.g., the equivalent 
diameter, D ), and the fluid properties of thermal conductivity, k, viscosity, p, specific 
heat, Cp, and density p 

The independent dimensional quantities to be used in the analysis are mass, M; 
length, L; time, O; and temperature, T. The variables, their symbols, and their 
dimensional equations are listed in Table 12. 



TABLE 12 

VARIABLES FOR CONVECTION 



VARIABLE 


SYMBOL 


DIMENSIONAL EQUATION 


Equivalent diameter 


D 


[L] 


Thermal conductivity 
of the fluid 


k 


[ML/0 3 T] 


Velocity of the fluid 


V 


[L/0] 


Density of the fluid 


P 


[M/L 3 ] 


Viscosity of the fluid 


P 


[M/L0] 


Specific heat at 
constant pressure 


c p 


[L 2 /0 2 T] 


Heat transfer 
coefficient 


h 


[M/0 3 T] 



There are seven physical quantities and four primary dimensions. We therefore expect 
that three dimensionless groups will be required to correlate experimental data. To 
determine these dimensionless groups, write as a product of each of the variables 
raised to an unknown power; 



50 



n = D a k b V c p d p e C p f h» 

and then substitute the dimensional formulas from the tabulation above: 

n = [L] a [ML/O 3 T] b [L/O] c [M/L 3 ] d [M/L0] e [L 2 /0 2 T] f [M/O 3 T]§ 

For n to be dimensionless, the exponents of each primary dimension must separately 
add to zero. Equating the sum of the exponents of each primary dimension to zero 



yields the following set of equations: 

b + d + e + g = 0 for M 

a + b + c — 3d — e + 2f = 0 forL 

— 3b — c — e — 2f — 3g = 0 for 0 

— b — f — g=0 for T 



Since there are seven unknowns but only four equations, three of the exponents can be 
arbitrarily selected for each dimensionless group. The only restriction is that each of 
the arbitrarily selected exponents be independent of the other. This requirement is met 
if the determinant formed by the coefficients of the remaining terms does not equal 
zero. 

Since the convection heat transfer coefficient is the dependent variable that we 
eventually want to evaluate, set its exponent g equal to unity in the First dimensionless - 
group, rij To simplify the algebraic manipulations, arbitrarily set c = d = 0. This 
means that in the First dimensionless group velocity and density will not appear. 

Solving the equations above simultaneously, we obtain a = 1, b = — 1, and e 
= F = 0. The First dimensionless group is therefore 

n 2 = hD/k 

which we recognize as the Nusselt number. 

For the second dimensionless group, we set g = 0 to avoid having the 
dependent variable h appear again. Arbitrarily, we let a = 1 and f = 0. Simultaneous 
solution of the equations above with these choices yields the second parameter, 

n 2 = VDp/p 

which we recognize as the Reynolds number, Re. 
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For the third dimensionless group, we wish to have the specific heat, which has 
not appeared in any previous group, to be included and therefore set f = 1. Similarly, 
we do not wish to have the heat transfer coefficient or the diameter appear again and 
therefore let a and g equal to 0. This yields the third dimensionless group, 

n 3 = CpM/k 

which we recognize as the Prandtl number, Pr. Thus the heat transfer coefficient in TT t 
can be functionally related to the dimensionless groups n 2 and Il 3 in the form 

Nu = f( Re, Pr ) 
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APPENDIX C 

UNCERTAINTY ANALYSIS 



The uncertainties for the major variable in the experiments were calculated in 
accordance with the method described by S. Kline and F. McClintoch (Ref. 9.). The 
estimates of the uncertainties were made conservatively. As a result, there is 
considerable confidence in the uncertainties as calculated. 

The following equations were used to calculate the uncertainties. 

1. BT b n b = [ (5T ,/T , ) 2 + (dT 2 /T 2 ) 2 ]'/ 2 

2. ac h /C h = [ (dm/m) 2 + (dC p /C p ) 2 ]‘/ 2 

3. dc = [ (3C h /C h ) 2 + (diyT,) 2 + ,(dT 2 /T 2 ) 2 + (dt ,/tj) 2 F 2 

4. dN tu /N tu = [ (dC r /C r ) 2 + (dm) 1 ] 1/2 

5. dU/U = [ (dN tu /N tu ) 2 + (dC h /C h ) 2 + (dA/A) 2 )>/ 2 

6. dV/V = [ (dm/m) 2 + (dp/p .) 2 + (dAf/A^) 2 J 1 / 2 • 

7. dh/h = [ (dU/U) 2 + (dr w /r w ) 2 )>/ 2 

8. dPr/Pr = ( (dC p /C p ) 2 + (dp/p) 2 + (dk/k) 2 ]'/ 2 

9. dRe/Re - ( (dp/p) 2 + (dV/V) 2 + (dD/D) 2 + (dp/p) 2 J 1 / 2 

10. dNu/Nu = [ (mdRe/Rc) 2 + (ndPr/Pr) 2 \' !l 

11. df/f = (mdRc/Rc) 

12. 3j/j = (m5Re/Re) 

The following quantities had their uncertainty calculated by dividing the estimated 
error in the quantity by the value of the quantity and these uncertainties are considered 
constant for the range of values in this experiment. 
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Quantity 



Uncertainty 



A 

A f 

D 

P 

k 

P 



m 

C„ 



YV 



0.003 

0 005 

0.006 

0.005 

0.003 

0.002 

0.008 

0.005 

0.003 



The following quantities had the worst case error 



Quantity 



Error 



T 

T 

t 



t 



2 



AT 



±0.2 

±0.2 

±0.2 

±0.2 

±0.3 



Then, the uncertainties for the following quantities were calculated using the equations 
above, 



Quantity 


Uncertainty 


T b 


0.008 


c h 


0.01 


c 


0.012 


N tu 


0.015 


u 


0.018 


h 


0.018 


V 


0.01 


Pr 


0.009 


Re 


0.013 


Nu 


0.009 
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f 

J 



0.005 

0.005 
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